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Abstract 

We report on a search for the parental molecular clouds of the TW Hya association (TWA), using CO 
emission and Nal absorption lines. TWA is the nearest young (~ 50 pe; ~ 10 Myr) stellar association, yet 
in spite of its youth, there are no detection of any associated natal molecular gas, as is the case for other 
typical young clusters. Using infrared maps as a guide, we conducted a CO cloud survey toward a region 
with a dust extinction of E{B — V)> 0.2 mag, or Ay > 0.6 mag. CO emission is detected toward three IR 
dust clouds, and we reject one cloud from the TWA, as no interstellar Na absorption was detected from 
the nearby Hipparcos stars, implying that it is too distant to be related. The other two clouds exhibit 
only faint and small-scale CO emission. Interstellar Nal absorptions of Hipparcos targets, HIP 57809, HIP 
64837, and HIP 64925 (at distances of 133, 81, and 101 pc, respectively) by these couds is also detected. 
We conclude that only a small fraction of the interstellar matter (ISM) toward the IR dust cloud is located 
at distance less than 100 pc, which may be all that is left of the remnant clouds of TWA; the remaining 
remnant cloud having dissipated in the last 1 Myr. Such a short dissipation timescale may be due to an 
external perturbation or kinematic segregation that has a large stellar proper motion relative to the natal 
cloud. 

1. Introduction of T Tauri stars at a range of distances to their nearest 

• ' molecular cloud (e.g., Neuhauser 1997). At this stage, 

'. Kastner et al. (1997) examined a few young stars near the process by which a T Tauri star becomes separated 
Q^the classical T Tauri star, TW Hya (which has become from its nursery molecular cloud is unclear; one plausible 
^^known as the TW Hya association) , and suggested that scenario postulates that the young stars are formed in rel- 
• 'they are nearby and of a young age. Later Hipparcos atively small clouds which are rapidly dissipated (Mizuno 
. J^piarallaxes measurements towards TW Hya and HD 98800 et al. 1998; Tachihara et al. 2001; Tachihara et al. 2005). 
^^yielded a distance of ~ 50 pc, and an age of ~ 10 Myr, In general however, young stars are found to be closely as- 
;_(a'lthough the cluster age and its dispersion are still unclear sociated with their natal molecular clouds. For example; 
.^('Tachihara et al. 2003). The stellar population number Taurus, Chamaeleon, and Lupus 3, whose ages are a few, 
of TWA has expanded after a ROSAT X-ray survey and 5, and a few xlO Myr, respectively, are all associated with 
other foUowup measurements (e.g., Webb et al. 1999), and a rich molecular cloud or envelope. An older example, the 
the current population of the association now comprises Pleiades (125 Myr) on the other hand, is not associated 
more than 20 stars (Zuckerman et al. 2001; Webb 2001). with any molecular gas. 

The proximity of TWA makes it an ideal source for A notable young exception is the ry Chamaeleontis clus- 
studying faint objects such as: extra-solar planets, proto- ter (~ 8 Myr) at ^ 97 pc, which has no surrounding 
planetary disks, and young brown dwarfs (Neuhauser et molecular gas despite its young age (Mamajek et al. 
al. 2000). Notably however, despite its young age, no 1999; Mizuno et al. 2001). Mamajek et al. (2000) proposed 
molecular cloud remnants from the formation of TWA that the rj Cha and TW Hya clusters were formed as parts 
have been found. of Sco-Cen OB association and later dispersed with rela- 

The ROSAT all sky survey and its follow-up optical ob- tively large proper motion. The proximity of TWA makes 
servations have been successful in finding a large number it better target for exploring this scenario in detail, as a 
- — r-— TT — - — search for remaining faint and small remnant clouds can 

* ESO run ID: 071.C-0568 1^ ^ -I 

t Present address: National Astronomical Observatory of Japan, °^ made more easily. 

2-21-2, Mitaka, Tokyo, 181-8588, Japan; k.tachihara@nao.ac.jp As the Star formation efficiency for typical Star forming 
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regions (SFRs) is much less than 50%, where the remain- 
ing molecular gas is returned to the ISM, the study of 
the dissipation of molecular gas is an important stepping- 
stone in the overall objective of understanding matter cy- 
cling in the Galaxy. The largest and most extensive molec- 
ular cloud survey in our Galaxy is the CfA CO survey 
(e.g.. Dame et al. 2001). However, this survey does not 
cover the distribution of the TWA members and it is not 
of sufficient sensitivity to detect the expectedly faint CO 
of the TWA molecular cloud remnant. 

2. Observations 

We have made a search for any remnant natal molec- 
ular cloud associated with TWA using the CO (J =1-0) 
emission line at 115.2712 GHz and Nal doublet absorp- 
tion lines at A = 5889.95 A and 5895.92 A. TWA subtends 
an approximate 40° x 30°, and we limit the survey size by 
using the the IRAS far-IR survey data as a guide map; 
Fig. 1 shows a dust extinction map of E{B — V) derived 
from IRAS 100 ^m and COBE/DIRBE 240 /xm, as demon- 
strated by Schlegel et al. (1998). We note that some in- 
frared cirrus clouds are located above the Galactic plane 
where TWA-member stars are distributed. However, these 
clouds do not appear to spatially correlate with the TWA- 
mcmber stars. The typical age of TWA stars is 20± 10 
Myr (Kastner et al. 1997; Jensen et al. 1998), and TWA 
members with proper motion of 1 km s~^ will travel ~ 10 
pc (~ 11 deg) since their formation, effectively separating 
them from the natal clouds. As such, we focus on rela- 
tively intense and large IR cirrus clouds, not biased onto 
spatial correlation of a TWA-member. We instead select- 
ing IR cirrus clouds whose E{B — V) peaks are greater 
than 0.2 mag (corresponding to Ay > 0.6 mag, assum- 
ing the empirical relation of Ay ~ Rv E{B ~ V) with 
Rv — 3.1). We identify one large and intense cloud, and 
two diffuse and extended clouds above b > 15 deg that con- 
tain several local peaks. Hereafter, we refer to these three 
prominent IR cirrus clouds shown in Fig. 1 as cloud A, B, 
and C. 

2.1. CO observations 

^^CO (J =1-0) observations were conducted with the 
4m NANTEN radio telescope at Las Campanas observa- 
tory in Chile. At 115 GHz, NANTEN has a beam FHWM 
of 2.6 arcmin, and a velocity resolution of 0.1 km s^^. 
Regions where the dust E{B — V) > 0.2 mag were quickly 
surveyed with a 4-arcmin grid spacing, to a rms sensitiv- 
ity of ^ 0.2 K. Resulting regions containing a possible CO 
detection were resurveyed to improve the sensitivity to 
< 0.1 K. These observations were made with a frequency- 
switching technique, using a throw interval of 13 - 20 
MHz depending on atmospheric condition. A position- 
switching technique was used instead for clouds with a 
Vlsr close to the telluric CO line. In total, 764 positions 
were observed. 
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2.2. Nal obsercations 

Wc check the physical association of the target molec- 
ular clouds with TWA through a distance comparison: 
TWA-member stars have a distance range of 30 to 120 pc 
(Frink 2001). Therefore, any background stars in the same 
line of sight of the clouds, and with a distance ofd> 120 pc 
should be measurably extincted by the candidate natal 
clouds. We selected 21 Hipparcos stars with a parallactic 
distances d < 150 pc, located in and around high dust col- 
umn density regions of cloud A-C. These were observed 
with the European Southern Observatory (ESO) 2.2m 
telescope using the FEROS (Fiber-fed Extended Range 
Spectrograph). The spectrum range of FEROS spans 3600 
to 9200 A, and has a resolving power or R = 48000 (2 pix- 
els resolution). The standard MIDAS pipeline was used 
for bias subtraction, flat-fielding, scattered light removal, 
sky background subtraction and wavelength calibration. 

3. Results 

3.1. CO survey 

The Integrated CO intensity map and the E{B — V) 
map for cloud A is shown in Fig 2, where we see relatively 
strong CO emission and its intensity distribution is gen- 
erally similar to that of dust with Ay > 0.84. Notably 
however, the strong dust peak at Z ^ 280.°6, 6 ~ 13.° 5 is 
weakly represented in the CO map. Conversely, the sec- 
ond strongest dust peas at I ^ 281.°5, b ^ 14.°5 is strongly 
represented in the CO map. In addition, we find another 
CO peak located between the two dust peaks. The peak 
CO integrated emission is 9.6 K km s~^, corresponding to 
a molecular hydrogen column density of 1.5 x 10^^ cm^^ 
(Hunter et al. 1997). This is roughly consistent with that 
derived from the dust Ay (~ 1.1 mag), assuming all the 
hydrogen is in the molecular form, i.e., iV(H2)/Ay = 10^^ 
cm~2 mag~^ (Bohlin et al. 1978). 

In contrast to cloud A, CO emission is detected from 
only small and weak regions in clouds B and C (Figs. 3 and 
4). In total, emission with W{CO) > 0.3 K km s~^ was 
measured at only 37 positions from total of 436 points, 
as shown in Figs. 3 and 4. Tabic 1 lists the properties 
of the brightest three W{CO) peaks each from clouds B 
and C. The dust maps for these two clouds shows a com- 
plex morphology, with several column density peaks hav- 
ing Ay > 0.6 mag. 

CO emissions are detected only toward the three of 
many peaks of cloud B. In contract, CO was detected to- 
wards all the three dust peaks of cloud C. The correlation 
between dust Ay and CO intensity is generally weak. The 
CO spectra towards both clouds B and C, shown in Figs. 3 
and 4, are generally complex and have a number of veloc- 
ity components. The systemic velocity of the CO emission 
components range between —9 to -t-1 km in cloud B, 
and from —8 to —5 km s~^ in cloud C. These velocities are 
slightly smaller than the mean Flsr of the TWA stars, at 
-t-2.8 km s-i (Torres et al. 2003). Additionally worthy of 
note, is that the Vlsr of the CO emission peaks are more 
dispersed than velocity dispersion of the TWA-member 
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Fig. 1. Far-IR dust Ay map derived from IRAS 100 ^m and COBE/DIRBE 240 ^m surveys (Schlegel et al. 1998). The TWA stars 
are shown by crosses with numbers. Three boxes designate the IR cirrus clouds named cloud A, B, and C from right to left. 



Table 1. Properties of CO peaks 



name 


Galactic 
longitude latitude 


R.A. 
(J2000.0) 


Decl. 
(J2000.0) 


W{CO) 
(K km s-i) 


peak Ti* 
(K) 


Vlsr* 
(km s-i) 


AV* 

(km s^"'^) 


N{R2y 

(1020 cm-2) 


Bl 


289.4667 


20.2667 


11 44 18.3 


-40 50 35 


0.32 


0.27 


-8.5 


0.9 


0.5 


B2 


291.4000 


18.1333 


11 50 49.8 


-43 22 38 


0.57 


0.32 


-0.8 


1.1 


0.9 


B3 


294.6000 


18.6667 


12 07 48.1 


-43 29 53 


0.93 


0.75 


1.0 


1.2 


1.5 


CI 


306.2667 


18.6000 


13 09 04.0 


-44 09 27 


2.02 


1.71 


-7.0 


0.9 


3.2 


C2 


306.8667 


19.0000 


13 12 02.3 


-43 42 51 


0.51 


0.26 


-5.0 


1.5 


0.8 


C3 


307.4667 


19.0667 


13 15 07.7 


-43 35 47 


0.90 


0.31 


-8.3 


2.4 


1.4 



"Derived by Htting to the gaussian proHles. 

tAssuming the conversion factor 7V(H2)/VK(CO) of 1.56 X 10^0 cm-^ (K km 3-^)^^ (Hunter et al. 1997). 



Molecular clouds in TW Hya 



3 




Fig. 3. Main: The dust E{B — V) map of cloud B. Contour intervals are 0.15+0.05 mag and the crosses indicated the position 
of the observed Hipparcos stars. Aside: Spectral maps for the three detected CO peaks. CO peak positions listed in Table 1 are 
labeled. The circle shows the position of HIP 57809 that shows interstellar Na absorption, while crosses are other HIP stars without 
Na detections. 



stars (1.9 km s^^ as a standard deviation). 

3.2. Optical spectroscopy 

To help distinguish the blended Nal stellar emission 
lines and interstellar absorption lines, we have defined an 
absorption detection according to the following three cri- 
teria: (a) Both of the Na doublet lines have two (partially) 
resolved components; (b) the interstellar lines are signifi- 
cantly narrower than the stellar ones; (c) the star is not a 
spectroscopic binary, according to many other absorption 
lines in the spectrum. Applying these criteria result in the 
detection of interstellar Na I D absorption of three targets; 
HIP 57809, HIP 64837, and HIP 64925 (Fig. 5). Targets 
HIP 57809 and HIP 64837 are towards clouds B and C, 
respectively, and HIP 64925 is located ^ 2 degrees away 
from cloud C. The parallactic distances for HIP 57809, 
HIP 64837 and HIP 64925 are 133, 81, and 101 pc, re- 
spectively. We do not detect absorption for any targets 



towards cloud A. 

The equivalent width, Vl^(Na) and Vlsr (calculated 
from the central wavelength) of the interstellar lines 
were estimated by fitting a double-component gaussian. 
To avoid confusion with atmospheric lines, we estimate 
VK(Na) from the D2 line only. After de-blended from the 
stellar absorption lines we estimate a W^(Na) of 0.01, 0.04, 
and 0.01 A and Vlsr of 13.6, -1.4, and -2.5 km s^^ for 
HIP 57809, HIP 64837, and HIP 64925, respectively. The 
optical spectral resolution relative to the radial velocity is 
6.1 km s~^, which is poorer than that obtained from the 
radio observations. 

4. Discussion 

4-.1. Distance to the clouds 

We use the Na absorption lines to constrain the distance 
to the three clouds. Nine Hipparcos stars with spectro- 
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Fig. 4. Main: The dust E{B — V) map of cloud C. Contour intervals are 0.15+0.05 mag and the crosses indicated the position 
of the observed Hipparcos stars. Aside: Spectral maps for the three detected CO peaks. CO peak positions listed in Table 1 are 
labeled. The circles show the positions of HIP 64837 and HIP 64925 that show interstellar Na absorption, while the crosse is HIP 
65178 without Na detection. 




5880 5885 5890 5895 
Wavelength (A) 



5900 5880 5885 5890 5895 5900 

Wavelength (A) 



5880 5885 5890 5895 5900 

Wavelength (A) 



5905 



Fig. 5. Optical spectra of HIP 57809, HIP 64837, and HIP 64925. The ordinate is the normalized relative intensity. They all show 
composite Na absorption lines that consist of wide stellar components and narrow interstellar ones. The intensities are in arbitrary 
scales. 
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Galactic longitude 

Fig. 2. Bottom: Dust E{B — V) map of cloud A. Contour 
intervals are 0.15+0.05 mag and the crosses indicate the po- 
sitions of the observed Hipparcos stars. Top: CO integrated 
intensity map of cloud A, over the region indicated in the 
dust map below. The solid lines indicate the limits of the CO 
surveyed region. Contour intervals are 1.7+1.7 K km s~^. 

scopic data are in the direction of cloud A, however none 
of these show evidence for absorption, thus the distance of 
the molecular cloud is likely to be greater than that of the 
nine stars, specifically. Da > 144 pc. The mean distance of 
TWA- members is much less than this, and so we conclude 
that cloud A is not associated with the TWA. The total 

cloud mass is derived is estimated by = 67.5 (^j^^-^^ ■ 

Determining the distances to clouds B and C is compli- 
cated by their fragmentary morphology. Four Hipparcos 
stars were observed towards cloud B and only one of 
them, HIP 57809 {d — 133 pc), shows interstellar absorp- 
tion. A nearby star {d = 92 pc), HIP 57860, does not 
show evidence for absorption. We can thus constrain 
the distance to cloud B (and more specifically the dust 
peak at I ~ 291.°8, b - 17.°7 east of the CO peak B2) to 
92 pc < DB2oast < 133 pc. The remaining two stars (HIP 
59120 and HIP 59257, at distances of 108 and 76 pc, re- 



spectively) do not show interstellar absorption, and so we 
can further constrain the distance estimate for the peak 
at / ~ 294.°5, 6 ~ 19 to Db3 > 108 pc. It is worth mention- 
ing that the velocity of the interstellar absorption toward 
HIP 57809 (Vlsr = +13.6 km s'^) is which is significantly 
larger than that of CO emission at the nearby peak B2 
(^lsr = —0.8 km s~^), suggesting two layers of ISM in 
the line of sight. Large difference in Vlsr of the three CO 
peaks (from —8.5 to +1.0 km s""'^) also imply multiple 
components. 

Cloud C shows similar level of complexity: Two 
Hipparcos stars are located close to this cloud (HIP 64837, 
d = 81 pc and HIP 65178, d== 128 pc) and only HIP 64837 
shows evidence for absorption. A third star, HIP 64925 
{d = 101 pc), also evidence for the absorption, however 
this star is somewhat more distance from the cloud (see 
Fig. 4). Together, this suggests that the cloud is at a dis- 
tance of Dc < 81 pc, although the nearby star HIP 65178 
at c? = 128 pc does not show absorption possibly due to 
fragmentary small scale cloud morphology. The veloci- 
ties of the Na absorption lines (Vlsr = — 1-4 km s~^) and 
CO emission lines (li,sR = —8.3 km s^^) are again sig- 
nificantly discrepant, although they are within the errors 
of the absorption resolution. We find, that the absorp- 
tion and emission spectra presented here suggest a com- 
plex multi-component spatial and kinematic distribution 
of clouds B and C. 

4-2. Natal clouds of the TW Hya association 

Considering the region encompassing TWA members 
that is not contaminated from the Galactic plane dust 
emission, 267° < I < 310° and 16° < & < 35°, we find that 
only ~ 3% or 23 deg^ shows a dust E{B-V)> 0.15 mag. 
The survey data presented in this report samples 1.94 deg^ 
(excluding cloud A), showing CO in emission at only 8.5% 
of the observed positions. We can extrapolate this detec- 
tion ratio to regions where E{B — V)> 0.15 mag, and es- 
timate that 1.9 deg^ will show CO in emission. Assuming 
a typical CO intensity of 0.5 K km s~^ and a molecular 
hydrogen column density of 8 x 10^^ cm~^, we futher esti- 
mate that the total molecular mass throughout the region 
defined above is ^ 3 Mq . 

We can also estimate the mass of original pre-stellar 
molecular cloud; assuming a total current stella mass of 
~ 16 Mq (Tachihara et al. 2003) and a star-formation ef- 
ficiency of 3% (Tachihara et al. 2002), yields a pre-stellar 
molecular cloud mass of 540 M©. As such, less than 1% 
of the original cloud mass remains since the most recent 
TWA star formation, approximately 1 Myr ago. This 
is a considerably short dissipation timescale, in compari- 
son to other nearby SFRs; Tachihara et al. (2001) noted 
that most T Tauri stars younger than 5 Myr found in 
Lupus and Chamaeleon, are within a few pc of the associ- 
ated molecular clouds. They also note that small isolated 
clouds may also have been dissipated rapidly by nearby 
OB associations. 

The present velocity dispersions of the candidate TWA 
natal molecular clouds are too high to be bound, suggest- 
ing that they may have been recently disturbed by some 
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external event. In the following, we present scenarios that 
may result in the observed and anomalously high velocity 
dispersions of the candidate TWA molecular clouds. 
4-2.1. Scenario I 

Careful analysis of kinematic data suggest that both 
the TWA and 77 Cha cluster formed from the Sco-Cen 
giant molecular cloud (GMC) at a distance of ~ 150 pc 
(Mamajek et al. 2000). Over the last 10 Myr, they 
have been moving away from the Sco-Cen complex and 
are presently separated from the Lower Centaurus Crux 
(LCC) subgroup by 50-100 pc. The Sco-Cen GMC was 
later evacuated by stellar winds and/or SN explosion 
forming the giant "Loop I bubble" Hi shell (de Geus 
1992). The progenitor of the 77 Cha cluster may have 
been dispersed by the same SN responsible for forming 
the Loop I bubble. 
4-2.2. Scenario II 

TWA is within 110 pc of the sun, and is inside a low- 
density plasma region that fills the local bubble, observed 
as isotropic soft X-ray radiation (e.g. Tanaka & Bleeker 
1977). The Pleiades moving group is a nearby OB asso- 
ciation that may have been the source of a SN explosion 
in the last ^ 1 Myr, or of 3 SNe explosions in the last 5 
Myr (Berghofer & Breitschwerdt 2002). The shock wave 
from these events may be responsible enhancing the dis- 
sipation of the TWA natal molecular clouds, however the 
Pleiades moving group is large, with a scale of ~ 100 pc, 
and as such the relative proximity of the SN events to 
TWA is unclear, furthermore, no candidate SNe of such 
an explosion has yet been identified. 

More accurate measurements of the proper motions, 
distances and radial velocities of both the TWA-member 
stars and candidate remnant molecular clouds, are cru- 
cial to refine and discriminate between these scenarios, as 
would search for any remnants of any SNs in the Pleiades 
moving group. 

5. Summary 

We have conducted a search for remnant clouds associ- 
ated with the formation of TW Hya association in both 
CO emission, and Na absorption. Wc find that CO emis- 
sion is detected towards three peaks of IR brightness, 
clouds A, B and C. Cloud A is relatively bright in CO 
emission and shows strong emission in IR, whereas clouds 
B and C are weaker by more than a factor of 3, and show 
a much more fragmentary distribution of CO emission. 
We reject the association of cloud A to the TWA, on the 
basis of a lack of a detection of Na absorption lines which 
indicate that it is much more distant than the TWA. Na 
absorption of Hipparcos targets help to constrain the dis- 
tances of clouds B and C as that only a small fraction of 
the ISM has a distance less than 100 pc. Using a statistical 
approach, we estimate the total amount of molecular gas 
in this region is less than a few solar masses, and conclude 
that most of the natal molecular clouds associated with 
TWA has already dissipated with an unusually fast dissi- 
pation time scale, enhanced perhaps, by an external and 
energetic disturbance, or kinetic segregation with proper 



motion of the member stars. 
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